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Under  this  DURIP  program,  we  designed,  procured,  and  integrated  a  LWIR  microgrid  polarimeter.  In 
addition,  we  performed  initial  training,  calibration,  and  verification  on  the  camera  system.  Multiple  intensity 
measurements  are  required  to  determine  the  polarization  state  of  light. ^  These  multiple  measurements  can  be 
made  sequentially  in  time  using  a  single  detector  (or  array  of  detectors  for  imaging  polarimetry),  or  they  can  be 
made  simultaneously  using  several  detector  arrays  located  at  different  physical  positions.  Numerous  strategies 
have  been  developed  to  affect  the  simultaneous  measurements,  as  discussed  in  a  recent  review  article.^ 

One  approach  that  has  been  developed  for  rapid  acquisition  of  polarimetric  data  is  to  place  polarization 
filters  directly  on  individual  detector  elements  of  the  FPA.^  This  advanced  FPA  technology  allows  simultaneous 
collection  of  data  sets  associated  with  each  of  the  polarizer  orientations  in  the  array.  The  specific  instrument 
that  we  are  working  with  is  a  640  x  480,  engineering-grade,  HgCdTe  FPA  manufactured  by  DRS  Sensors  and 
Targeting  Systems  operating  in  the  band  8-10  fim.  The  FPA  has  an  integrated  wire  grid  polarizer  array  that 
is  aligned  and  bonded  to  the  FPA.^  The  polarizer  arrangement  forms  an  alternating  pattern  of  0,  45,  90,  and 
135  degree  orientations  as  shown  in  Fig.  1,  and  the  measured  intensity  data  is  processed  to  estimate  the  first 
three  (sq,  si,  S2)  Stokes  parameters.^  In  principle  only  three  measurements  are  necessary,  but  the  fourth 
measurement  provides  robustness,  reduces  noise,  aids  in  eliminating  the  effects  of  dead  pixels  in  the  FPA^  and 
fits  well  with  the  rectangular  nature  of  the  FPA  pixel  pattern. The  polarization  estimate  is  reconstructed 
at  the  intersection  of  the  four  pixels  that  compose  a  single  sub-quad  of  pixels  or  a  “super  pixel”  as  depicted 
in  Fig.  2  using  one  of  several  interpolation  strategies. These  polarized  FPAs  allow  for  direct  replacement  of 
existing  thermal  FPAs,  augmenting  the  thermal  imagery  by  providing  additional  information  in  real  time. 

The  FPA  was  purchased  from  DRS  Sensors  and  Targeting  Systems.  The  cooler  was  purchased  from  L3 
Technologies  (Cincinnati  Electronics).  The  system  was  integrated  by  SE-IR  in  Santa  Barbara,  CA.  We  are 
using  a  Janos  E/2,  50-mm- focal  length  LWIR  lens,  and  we  have  integrated  a  band-pass  filter  with  a  pass  band 
of  7.9  -  9.9  microns. 


O  Reconstruction  Point  {m  \  n’) 
Microgrid  Superpixel 


Figure  1.  Division  of  Focal  Plane  (DoFP)  polarimeters  operate  by  placing  a  micropolarizer  in  front  of  each  pixel  that 
analyzes  a  different  polarization  state.  Stokes  vectors  are  estimated  at  interpolation  points  as  shown  in  the  figure.^ 


Actual  Microgrid  Pixel 


Figure  2.  Super-pixel  of  the  microgrid  polarimeter. 


Figure  3.  LWIR  microgrid  camera  (center)  with  the  control  electronics  (right)  being  calibrated  with  a  blackbody  (blue 
box  on  left). 


the  SE-IR  system  came  with  their  CAM-IRA  software  package.  We  use  Cam-IRA  for  diagnostics  and  testing, 
but  we  have  also  written  our  own  custom  libraries  that  we  use  for  experimentation  and  data  collections.  The 
camera  system  was  integrated  onto  a  portable  workstation  so  it  could  be  taken  out  of  the  lab  for  testing.  The 
camera  is  shown  in  our  laboratory  in  Fig.  3. 

Example  polarimetric  imagery  taken  both  inside  the  lab  and  outside  the  lab  are  shown  in  Fig.  4.  A  recently 
published  paper  from  Optics  Letters  it  appended  to  this  report  as  Appendix  A  that  includes  data  taken  with 
the  DURIP  polarimeter. 


Figure  4.  Left:  Image  of  a  warm  emissive  sphere  taken  in  our  laboratory  at  the  College  of  Optical  Sciences.  Right; 
Image  of  the  Flandrau  Science  Center  on  the  campus  of  the  University  of  Arizona.  Taken  from  the  roof  of  the  Meinel 
Building,  College  of  Optical  Sciences. 
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Microgrid  polarimeters  operate  by  integrating  a  focal  plane  array  with  an  array  of  micropolarizers.  The 
Stokes  parameters  are  estimated  by  comparing  polarization  measurements  from  pixels  in  a  neighborhood 
around  the  point  of  interest.  The  main  drawback  is  that  the  measurements  used  to  estimate  the  Stokes  vec¬ 
tor  are  made  at  different  locations,  leading  to  a  false  polarization  signature  owing  to  instantaneous  field-of- 
view  (IFOV)  errors.  We  demonstrate  for  the  first  time,  to  our  knowledge,  that  spatially  band  limited  polar¬ 
ization  images  can  be  ideally  reconstructed  with  no  IFOV  error  by  using  a  linear  system  framework.  ©  2009 
Optical  Society  of  America 
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Recent  years  have  seen  the  emergence  of  passive  po- 
larimetric  imagery  as  an  important  discriminating 
tool  in  many  areas  of  sensing.  Polarization  has  been 
used  to  help  in  target  detection,  defeat  clutter,  and 
see  through  mildly  turbid  media,  among  other  appli¬ 
cations.  A  recent  review  paper  [1]  contains  an  in- 
depth  discussion  of  many  of  these  applications. 

Polarized  incoherent  light  is  often  described  in 
terms  of  the  Stokes  vector, 

S  = [Sq  Si  S2  §3]^ 

=  Uh-^Iv  Ir-Iv  ^45-^35  (1) 

where  Ijj,  ^45?  ^135?  Ip  are  the  intensities 

observed  through  horizontal,  vertical,  45°,  135°,  left- 
circular,  and  right-circular  polarizers,  respectively. 
Microgrid  polarimeters  combine  an  array  of  micropo¬ 
larizers  with  a  focal  plane  array  (FPA)  to  make  a  spa¬ 
tially  modulated  measurement  of  the  polarized  inten¬ 
sity  [2] .  The  most  common  configuration  is  as  shown 
in  Fig.  1,  where  four  linear  polarization  states  are  in¬ 
terlaced  to  estimate  the  three  linear  polarization 
Stokes  parameters  Sq,  Si,  and  S2.  However,  systems 
have  been  proposed  where  other  polarization  mea¬ 
surements  are  made,  including  full  polarimeters  that 
also  measure  the  circularly  polarized  component 
[3,4]. 

Microgrid  polarimeters  are  attractive  because  they 
are  rugged  and  inherently  optomechanically  aligned. 
The  largest  drawback  of  these  systems  is  that  the 
Stokes  parameters  are  estimated  using  polarization 
measurements  from  different  locations  in  the  FPA. 
When  there  are  intensity  variations  in  addition  to  po¬ 
larization  variations  across  the  scene,  the  different 
instantaneous  fields  of  view  (IFOVs)  of  neighboring 
pixels  complicate  the  reconstruction  process  and  are 
the  primary  source  of  false  polarization  signatures 
[5]. 

All  studies  known  to  us  reconstruct  the  polariza¬ 
tion  estimates  from  the  raw  microgrid  image  by  us¬ 
ing  spatially  local  interpolation  kernels  [1,2,4,6,7].  A 


recent  study  demonstrated  that  the  false  polarization 
signatures  in  these  cases  are  because  of  an  incom¬ 
plete  demodulation  of  the  intensity  image  [5] . 

Microgrid  polarimeters  are  a  class  of  “snapshot” 
polarimeters  that  determine  all  of  the  desired  polar¬ 
ization  information  from  a  single  measurement  of  in¬ 
tensity  on  a  single  FPA.  Prismatic  polarimeters 
[8-10]  are  another  class  of  snapshot  polarimeter  that 
use  two  spatially  varying  multiorder  wave  plates  (bi- 
refringent  prisms)  to  introduce  polarization- 
dependent  carrier  fringes  that  are  modulated  by  the 
image  intensity,  thereby  creating  polarization- 
dependent  sidebands  in  the  spatial  frequency  plane. 
Comparing  the  results  of  Ratliff  et  al.  [5]  for  micro¬ 
grid  instruments  to  those  from  prismatic  polarim¬ 
eters  led  us  to  the  present  Letter,  which  uses  a  linear 
system  framework  to  create  ideal  microgrid  interpo¬ 
lation  kernels  that  can  perfectly  reconstruct  band 
limited  polarization  imagery 

We  develop  the  theory  here  for  the  specific  configu¬ 
ration  shown  in  Fig.  1;  however,  the  theory  is  easy  to 
extend  to  an  arbitrary  micropolarizer  layout.  We  will 
assume  that  the  pixels  in  the  FPA  are  ideal  point 
samplers;  the  effects  of  finite  pixel  size  and  shape  can 
be  included  using  the  standard  sampling  theory  The 
image  is  described  by  the  spatially  varying  Stokes 


m 


O  Reconstruction  Point  {m  \  n  ’) 
Microgrid  Superpixel 


Fig.  1.  (Color  online)  Typical  microgrid  polarimeter  uses 
four  interlaced  linear  polarizer  orientations  to  estimate  Sq, 
Si,  and  S2. 
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vector  image  Si(m,n).  The  intensity  measured  by  the 
(m,n)th  pixel  can  be  written  as 

Iim,n)  =  SAim,nf  ■Si(m,n),  (2) 


where  is  the  analyzer  Stokes  vector  of  the 

(m,n)th  pixel  that  is  expressed  as 


:  {m,n) 


|(cos(m7r)  +  cos(7i7r)) 
|(cos(m7r)  -  cos(7i7r)) 
0 


(3) 


Substituting  Eq.  (3)  into  Eq.  (2)  yields 

I{m,n)  =  SQ{m,n)  +  |cos(m7r)[si(7n,n)  +  S2(m,n)] 

+  |cos(7i7r)[si(m,n)  -  S2(m,n)].  (4) 

Next  we  take  the  discrete-space  Fourier  transform 
(DSFT)  of I(m,n).  We  assume  that  the  Stokes  param¬ 
eter  images  are  spatially  band  limited  enough  to 
avoid  aliasing  (this  condition  will  be  made  explicit  be¬ 
low),  and  we  consider  only  the  positive  frequencies  as 

I(i,  7})  =  v)  +  v)  +  S2(^-  77)] 

+  I) -S2(^,77- |)]>  (5) 

where  ^  and  rj  are  the  horizontal  and  the  vertical  spa¬ 
tial  frequencies  (in  cycles  per  pixel),  respectively  In 
Eq.  (5)  there  is  a  baseband  signal  carrying  Sq,  a  hori¬ 
zontal  sideband  signal  at  that  represents  Si-\-S2, 
and  a  vertical  sideband  at  77=  |  that  represents  Si 
-^2. 

To  demonstrate  an  ideal  reconstruction  we  use  the 
simple  polarized  scene  described  by 

(6) 


Si{n,m)  =  ^SQ(n,m)eri(2m),  (7) 


Fig.  2.  (Color  online)  Fourier  transform  of  I{m,n)  for  the 
excitation  defined  by  Eqs.  (6)-(8).  The  outer  shaded  regions 
represent  the  high-pass  filters  used  to  reconstruct  Si  and 
S2.  The  central  shaded  area  is  the  low-pass  filter  used  to  re¬ 
construct  Sq. 


pixels  that  contact  that  node  [5] .  The  NLPN  interpo¬ 
lation  is  the  most  widely  used  reconstruction  method 
for  microgrid  polarimeters  [2,4,6,7]. 

The  general  notion  in  the  polarization  imaging 
community  is  that  the  NLPN  interpolation  can  be 
used  to  create  good  reconstructions  by  simply  choos¬ 
ing  the  point  spread  function  of  the  optics  to  band 
limit  the  image.  A  second  notion  is  that  microgrid  po¬ 
larimeters  can  only  reconstruct  accurately  scenes 
that  have  features  that  are  constant  across  a  single 
4x4  superpixel.  Figure  3  clearly  demonstrates  that 
these  notions  are  false.  What  error  is  present  in  Figs. 
3A  and  3B  is  because  the  images  defined  in  Eqs. 
(6)-(8)  are  not  in  fact  band  limited,  resulting  in  small 
aliasing  effects.  Other  more  complicated  interpola- 


S2{n,m)  =  ^SQ(n,m)eri(2n) .  (8) 

The  Fourier  transform  of I(m,n)  for  this  excitation  is 
shown  in  Fig.  2.  We  see  both  the  baseband  signal  that 
corresponds  to  the  Fourier  transform  of  Sq  and  the 
horizontal  and  the  vertical  sidebands  predicted  in 
Eq.  (5). 

We  can  use  Fig.  2  to  form  a  quantitative  descrip¬ 
tion  of  the  band  limit  requirements.  There  are  many 
ways  that  the  spatial  frequency  information  among 
Sq,  Si,  and  S2  can  be  distributed  to  allow  a  perfect  re¬ 
construction.  However,  if  we  assume  that  Sq  is  spa¬ 
tially  band  limited  to  Wq  and  Si  and  S2  are 

spatially  band  limited  to  a  sufficient 

condition  to  avoid  aliasing  is  that  Wq-\-Wi<^. 

The  reconstruction  errors  for  Sq  and  Si  using  the  fil¬ 
ters  indicated  in  Fig.  2  are  shown  in  Fig.  3  for  our 
method  (Fig.  3A)  and  the  standard  nearest  like- 
polarization  neighbor  (NLPN)  interpolation  strategy 
(Fig.  3B).  The  NLPN  method  reconstructs  the  Stokes 
vector  estimates  at  each  node  in  Fig.  1  using  the  four 


-3  -3 

x10  x10 


Fig.  3.  Sq  (left  column)  and  Si  (right  column)  error  distri¬ 
butions  for  the  ideal  reconstruction  strategy  presented  here 
(top  row)  and  NLPN  interpolation  (bottom  row). 
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tors  than  the  NLPN  method  can  do  a  better  job  of  re¬ 
constructing  the  polarization  information  than  the 
NLPN  [5],  but  the  improvement  in  those  cases  is  by 
less  than  1  order  of  magnitude  (data  not  shown).  This 
occurs  because  the  interpolation  kernel  for  the  NLPN 
is  a  2  X  2  rectangular  window.  This  interpolation  ker¬ 
nel  will  have  a  corresponding  sinc-function  filter  in 
the  frequency  domain  that  allows  significant  spectral 
leakage  of  the  sidelobe  signals  into  the  reconstructed 
baseband  signal,  and  vice  versa.  This  can  be  inter¬ 
preted  as  a  form  of  polarimetric  aliasing. 

A  comparison  between  the  NLPN  and  the  linear 
system  method  presented  here  is  shown  for  real  data 
in  Fig.  4.  The  object  here  is  a  45  cm  diameter  brass 
sphere  painted  with  krylon  flat  black  paint  and 
heated  with  an  incandescent  light  bulb  to  50 °C.  The 
sphere  is  located  in  a  laboratory  with  a  nominal  room 
temperature  of  23°C.  It  is  expected  that  the  predomi¬ 
nant  polarization  signature  around  the  sphere  will  be 
p  polarized,  since  the  emitted  radiation  dominates 
the  image.  The  degree  of  polarization  is  expected  to 
be  maximal  near  the  edges  where  the  photons  leave 
the  sphere  near  grazing,  and  the  degree  of  polariza¬ 
tion  near  the  center  of  the  sphere  should  be  approxi¬ 
mately  zero.  Figures  4C  and  4D  are  colored  using  a 
transform  that  maps  the  angle  of  polarization  into 
color  hue  and  the  degree  of  linear  polarization  into 
color  saturation  [11].  To  mitigate  the  effects  of  alias¬ 
ing  we  used  Hamming  windows  for  our  frequency  do¬ 
main  filters. 

Microgrid  polarimetric  imagers  were  an  extremely 
popular  area  of  research  owing  to  their  ruggedness 


Fig.  4.  (Color  online)  Long-wave  IR  (8-10  pm)  polarimet- 
ric  imagery  taken  using  a  microgrid  polarimeter  described 
elsewhere  [12].  The  target  is  a  spherical  gray  body  at  50 °C. 
A,  Calibrated  microgrid  image  with  modulation  present.  B, 
Fourier  transform  of  I{m,n).  C,  Reconstruction  using  Eq. 
(5).  D,  Reconstruction  using  NLPN  interpolation  [5]. 


and  inherent  optomechanical  alignment,  but  they 
lost  favor  because  of  the  polarimetric  aliasing  prob¬ 
lems  near  contrast  edges  in  images  apparent  in  Fig. 
4D.  In  this  Letter  we  have  presented  a  reconstruction 
paradigm  that  eliminates  these  artifacts,  making  mi¬ 
crogrid  instruments  competitive  with  other  polari¬ 
metric  modalities  [1].  By  operating  in  the  frequency 
domain,  it  is  possible  to  perfectly  reconstruct  spa¬ 
tially  band  limited  polarization  scenes  with  no  IFOV 
error  for  the  first  time  known  to  the  authors.  In  cer¬ 
tain  instances  it  might  be  preferable  to  perform  the 
reconstruction  without  taking  a  Fourier  transform. 
For  example,  several  researchers  have  developed 
methods  to  perform  the  polarimetric  reconstruction 
in  hardware  at  the  FPA  by  making  local  arithmetic 
computations  [2,6].  The  results  derived  here  can  be 
used  to  derive  spatially  local  interpolation  kernels 
that  have  greatly  improved  error  performance  over 
the  methods  currently  employed  [5] .  Finally,  going  to 
the  frequency  domain  allows  us  to  consider  the  appli¬ 
cation  of  superresolution  methods  that  can  exploit 
aliasing  to  improve  the  resolution  beyond  the  super¬ 
pixel  or  even  individual  pixel  levels.  This  is  an  active 
area  of  current  work  [13]. 
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